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SUMMARY 

The existing conventional methods of correlating and predicting friction coefficients 
for laminar and turbulent flow, where the physical properties and density do not vary 
greatly, are shown to give friction coefficients that are in poor agreement with the mea- 
sured values when there are large variations in the physical properties, that is, large ra- 
tios of surface to fluid-bulk temperature. 

The local and average friction coefficients used were measured by seven investiga- 
tors for laminar and turbulent flow of helium, hydrogen, nitrogen, carbon dioxide, and 
air through smooth tubes. Inside diameters varied from 0. 115 to 0. 569 inches (0.292 to 
1. 445 cm); ratios of surface to fluid-bulk temperature ranged from 0. 35 to 7. 35; and 
modified surface Reynolds numbers ranged from 170 to 550 000. These data were used to 
determine the best methods of correlating and predicting local friction coefficients for ra- 
tios of distance from entrance of test section to inside diameter of test section (x/D) from 
16 to 113 and average friction coefficients for ratios of length to diameter (L/D) from 
21 to 200. 

The recommended correlation equation for modified surface Reynolds numbers less 

than 3000 is f/2 = 8/Re„, where f/2 is half friction coefficient and Re_ is the modified 

s s 

surface Reynolds number. For modified surface Reynolds numbers of 3000 or greater, 
the recommended correlation is f/2 = (0.0007 + 0. 0625/Reg’ (Tjj/Tg)®*®, where T^ 
and T are the bulk and surface temperatures, respectively. The foregoing smooth tube 
relations also correlated laminar and turbulent friction coefficients for flow between par- 
allel plates. 



INTRODUCTION 


There has long been a need for a means of correlating both laminar and turbulent 
friction coefficients for gases, with large variations in the physical properties, flowing 
through smooth tubes. Probably the most widely used method of correlating and predict- 
ing friction coefficients for turbulent flow is that of Karman-Nikuradse, in which the fric- 
tion coefficient and Reynolds number are evaluated at the film temperature Tj, as postu- 
lated by Humble, Lowdermilk, and Desmon (ref. 1): 



All symbols are defined in appendix A 

This method worked well for Reynolds numbers greater than 20 000, for surface to 
fluid-bulk temperature ratios T g /T b less than 2. 5, and for a length to diameter ratio of 
60. Recent experiments present friction coefficients that cannot be correlated by equa- 
tion (1). In these experiments, Taylor (refs. 2 and 3) measured average friction coeffi- 
cients for helium and hydrogen with wall temperatures up to 5600° R (3110° K) and T g /T b 
from 1.0 to 4. 1, and Perkins and Worsoe-Schmidt (ref. 4) measured local friction coef- 
ficients for nitrogen with T g /T b from 1.4 to 7. 35. Figure 1 shows that the measured 
friction coefficient can be as much as three times the value predicted by equation (1). 
Perkins, et al. found that some of their average friction data (T g /T b < 2. 46) could be 
correlated with T g /T b raised to the 0. 6 power. 

Dalle-Donne (ref. 5) stated that the bulk friction coefficient can be predicted by the 

use of the modified surface Reynolds number, for both laminar and turbulent flow, for 

1 < T /T, <2.2. When the correlation method of reference 5 was extended to T /T. 

s b — ® ° 

as great as 7.35, the experimental data fell as much as 70 percent below the Karman- 

Nikuradse line, as shown in figure 2. 

Other experimental data are now available. Davenport (ref. 6) presented heating data 
for helium and nitrogen in the laminar region. For turbulent flow, McEligot (ref. 7) and 
Magee (ref. 8) give heating data for helium and air, Wolf (ref. 9) gives both heating and 
cooling data for air and carbon dioxide, and McCarthy and Wolf (ref. 10) give heating data 
for helium and hydrogen. The data of reference 10 were not used in the present investi- 
gation because the magnitude of the momentum pressure drops approached (and sometimes 
exceeded) that of the total pressure drop (see appendix B) . The range of conditions cov- 
ered by the various investigators is shown in table I. 
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Half friction coefficient, f f /2 


.1 



Modified film Reynolds number Re f = (GD/pj)(T b /T f ) 


Figure 1. - Variation of local and average friction coefficients with modified film Reynolds number. Viscosity, and density evaluated at film temperature. 
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Modified surface Reynolds number, Re s = (GD/(j $ )(T b /T s ) 

Figure 2. - Variation of local and average friction coefficients with modified surface Reynolds number. Density in friction coefficient evaluated at bulk 
temperature; viscosity in Reynolds number evaluated at surface temperature. 



Presented herein are methods of correlating both local and average friction coeffi- 
cients for laminar and turbulent flow of gases through smooth tubes with ratios of surface 
to fluid-bulk temperature from 0. 35 to 7. 35. These friction coefficients were obtained 
for helium, hydrogen, nitrogen, air, and carbon dioxide and represent the work of seven 
investigators . 

The equations used in this investigation to correlate the friction data can easily be 
used to predict friction coefficients for designers working within the range of test condi- 
tions presented in table I. 


METHOD OF CORRELATION 

Dalle-Donne was able to correlate friction coefficients for both laminar and turbulent 
flow of air and helium with small variations in the physical properties (1 < T g /T b < 2.2) 
by using the bulk friction coefficient and a modified surface Reynolds number. In this 
report, their correlation method is compared with data for extreme variations in physical 
properties due to either heating (T g /T^ > 1) or cooling (Tg/T^ < 1). 

The equations presented herein were also compared with the data for laminar and 
turbulent friction coefficients for flow between parallel plates reported by Slaby, Maag, 
and Siegel (ref. 11). 

The relation recommended by Dalle-Donne for laminar flow with 1 < T /T, <2.2 

s b 

is 


2 Re g 

Equation (2) was used in this investigation to correlate friction data for helium, hydrogen, 
nitrogen, and air with 1 < T g /T b <4.1 with good results for modified surface Reynolds 
numbers up to 3000 (bulk Reynolds numbers up to 22 000), as shown in figure 2. 

Dalle-Donne also found that the use of the bulk friction coefficient and the modified 
surface Reynolds number correlated their friction data for turbulent flow with 1 < Tg/T^ 
<2.2 to the Karman-Nikuradse line with good results. In the present investigation, the 
Koo, Drew, and McAdams relation (ref. 12) was used because of its simplicity and its 
close agreement with the Karman-Nikuradse relation. The Koo, Drew, and McAdams 
relation for bulk Reynolds numbers from 3000 to 3 000 000 is given as 


-= 0.0007 + 
2 


0.0625 


Re 


0.32 


( 3 ) 
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In the present investigation, a modified surface Reynolds number was used in place of the 
bulk Reynolds number and gives the following relation: 


-= 0.0007 + 
2 


0.0625 


Re 


0.32 


(4) 


When equation (4) was used to predict friction coefficients for 0. 35 < T g /T b <7.35 
and for modified surface Reynolds numbers from 3000 to 550 000, the measured values 
deviated from the predicted values by as much as -70 percent for the runs with heat addi- 
tion (1. 0 < T g /T b < 7. 35) and by as much as 70 percent for the runs with heat removal 
(0. 35 < T g /T b < 1. 0), as shown in figure 2. Since the deviation appeared to be a function 
of T s /T b , the ratio of the friction coefficient calculated by equation (4) to the measured 
friction coefficient was plotted as a function of T /T fe (fig. 3). The slope of the line 
drawn through the data points is 0. 5 and passes through (f/2) . /(i/2) = 1 at 

T g /T b =1, as it should since equation (3) is equal to equation (4) at a T g /T b = 1. Fig- 
ure 3 indicates that equation (4) should be modified with the square root of T g /T b as 
follows : 


f 

2 


0.0007 


\ /rp \0. 5 
[ 0. 0625\ / p\ 

r< 32 /VV 


(5) 



r 


Figure 3. - Variation of ratio of friction coefficient calculated by equation (4) to 
measured local and average friction coefficient with ratio of surface to fluid-bulk 
temperature for modified surface Reynolds numbers of 3000 and greater. 
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Equation (5) predicts friction coefficients within ±10 percent for both heating and cool- 
ing of gases with modified surface Reynolds numbers from 3000 to 550 000 (bulk Reynolds 
numbers from about 5400 to 187 000, the limit of the experimental data). 

Since figure 3 shows an effect of Tg/T^ on the friction coefficient for modified sur- 
face Reynolds numbers of 3000 and greater, a similar plot was made for modified surface 
Reynolds numbers less than 3000. Figure 4 shows the ratio of the friction coefficient cal- 
culated by equation (2) to the measured friction coefficient as a function of Tg/T^ for a 
modified surface Reynolds number less than 3000. It appears from, figure 4 that there is 
no discernable effect of Tg/T^ on the friction coefficient for a modified surface Reynolds 
number less than 3000. 
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Figure 4. - Variation of ratio of friction coeffi- 
cient calculated by equation (2) to measured 
local and average friction coefficient with 
ratio of surface to fluid-bulk temperature for 
modified surface Reynolds numbers less than 
3000. 


DISCUSSION OF RESULTS 


The relations shown in equations (2) and (5) correlate local friction coefficients for 
x/D from 16 to 113 and average friction coefficients for L/D from 21 to 200 (with 
0. 35 < Tg/T^ <7.35 and modified surface Reynolds numbers from 170 to 550 000) for 
both vertical and horizontal smooth tubes. The tube inside diameters range from 0.115 
to 0. 569 inch (0. 292 to 1. 445 cm). 

In figure 5 the experimental friction coefficients are shown as a function of the modi- 
fied surface Reynolds number for low Reynolds number flow (Re„ < 3000) . Of the 109 ex- 

s 

perimental points shown, 84 percent of them fall within ±20 percent of the correlation line 
represented by equation (2) . 
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Figure 5. - Correlation of local and average friction coefficients 
for modified surface Reynolds number less than 3000. Den- 
sity in friction coefficient evaluated at bulk temperature; vis- 
cosity in Reynolds number evaluated at surface temperature. 
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Figure 6. - Correlation of local and average friction coefficients for modified surface Reynolds 
numbers of 3000 and greater. Density in friction coefficients evaluated at bulk temperature; 
viscosity in Reynolds number evaluated at surface temperature. 
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Figure 6 shows the experimental friction coefficients as a function of the modified 
surface Reynolds number for turbulent flow (Re g of 3000 and greater). There are 423 ex- 
perimental data points for both heating and cooling of the gas; 97 percent of the data 
points are within ±15 percent, and 90 percent fall within ±10 percent of the correlation 
line defined by equation (5). 

The data in figure 5 indicate a lack of a transition region, except for the few points 
with a Tg/T^ <1.3 that fall below the correlation line. Thus, it appears that the large 
variation in the physical properties and density tends to make the flow turbulent at a modi- 
v fied surface Reynolds number of 3000 or more. 

Slaby, Maag, and Siegel (ref. 11) measured laminar and turbulent friction coeffi- 
* cients for flow between parallel plates and attempted to correlate the data by the generally 

accepted method of evaluating physical properties and density at a film temperature. The 
results are shown in figure 7. The agreement between measured and predicted friction 
coefficients is poor for modified film Reynolds numbers below 3500. The agreement is 
better at higher Reynolds numbers. Better agreement between measured and predicted 
friction coefficients occurs when equations (2) and (5) are used. Figure 8 shows the cor- 
relation of measured parallel-plate friction coefficients using equations (2) and (5) . 


20xl0 -3 



Figure 7. - Correlation of average friction coefficients for 
laminar and turbulent flow between parallel plates with 
physical properties and density evaluated at film tempera- 
ture. Experimental data from reference 11. 
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Figure 8. - Correlation of average friction coefficients for laminar and 
turbulent flow between parallel plates with use of equations (2) and 
(5). Experimental data from reference 11. 


CONCLUSIONS 

For ratios of surface to fluid-bulk temperature of 0. 35 to 7. 35, modified surface 
Reynolds numbers from 170 to 550 000, and ratios of length to diameter from 21 to 200, 
the following conclusions may be made for flow through smooth circular tubes. These 
conclusions also appear to apply to flow between parallel plates. 

1. For modified surface Reynolds numbers below 3000, both local and average fric- 
tion coefficients for ratios of surface to fluid-bulk temperature of 1. 0 to 4. 1 can be pre- 
dicted within ±20 percent by the Dalle-Donne relation f/2 = 8/Re g where f/2 is the 
half friction coefficient and Re c is the modified surface Reynolds number. 

2. For modified surface Reynolds numbers of 3000 and greater, the flow appears 

to be fully turbulent. Both local and average friction coefficients for ratios of surface 
to fluid-bulk temperature from 0. 35 to 7. 35 can be predicted within ±10 percent by the 
relation f/2 = (0.0007 + 0. 0625/Re g ‘ (T^/Tg) 0, 5 , where f/2 is the half friction coef- 

ficient, Re g is the modified surface Reynolds number, and and T g are the bulk and 
surface temperatures, respectively. 

3. For ratios of surface to fluid-bulk temperature greater than approximately 1.3, 
there appears to be no transition region for friction coefficients. 

4. The modified surface Reynolds number determines if the laminar flow equation or 
the turbulent flow equation should be used to predict friction coefficients. 
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5. For large ratios of surface to fluid-bulk temperature, the use of the reference- 
temperature concept (evaluating physical properties and density in both the friction coef- 
ficient and Reynolds number at some reference temperature [T^ s X(T g - T^) + T^] is not 
applicable to friction coefficients for flow through smooth tubes or flow between parallel 
plates. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 25, 1967, 

122-28-02-16-22. 
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APPENDIX A 


SYMBOLS 

D inside diameter of test section, ft; m 

o 

dp/dx pressure gradient in axial direction, (N/m )/m 
f/2 half friction coefficient, S C P^ V APf r /4(L/D)G^ 
fj/2 half film friction coefficient, (f/2)(Tj/T^) 

G mass flow per unit cross-sectional area, lb/(hr)(ft ); kg/(sec)(m ) 

g c conversion factor, 4. 17x10^ (lb mass) (ft)/ (lb force) (sec'*); 9. 8 (kg)(m)/(N)(sec^) 

L test section length, ft; m 

p static pressure, lb/ft ; N/m 

Ap static pressure drop across test section, lb/ft ; N/m 

R specific gas constant 

Re^ bulk Reynolds number, GD/p^ 

Re^ modified film Reynolds number, (GD/ p^) (T^/T^) 

Re_ modified surface Reynolds number, (GD/p UT^/T J 
bulk temperature, °R; °K 
Tj film temperature, (T g + T^J/2, °R; °K 
T surface temperature, °R, °K 

reference temperature, °R; °K 
t static temperature, °R; °K 

V velocity, ft/hr; m/sec 

X parameter used in reference temperature equation 

x distance from entrance of test section, ft; m 

o 

p^ absolute viscosity of gas at T^, lb/(hr)(ft); (N)(sec)/m 

2 

p f absolute viscosity of gas at T f , lb/(hr)(ft); (N)(sec)/m 

2 

p o absolute viscosity of gas at T , lb/(hr)(ft); (N)(sec)/m 

s s 

p olr average density of gas (eq. (B5)) 

av 

p bay average bulk density of gas based on average static pressure and average static 
bulk temperature (eq. (B4)), lb/ft kg/m^ 
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2 2 

t shear stress at wall, lb/ft ; N/m 

w 

Subscripts: 

calc calculated value 

exp measured value 

fr friction 

mom momentum 

1 inlet 

2 outlet 


APPENDIX B 


CALCULATION OF FRICTION COEFFICIENTS 
Average Friction Coefficients 

The experimentally determined average friction coefficients for smooth tubes pre- 
sented by Taylor (refs. 2 and 3), Wolf (ref. 9), McCarthy and Wolf (ref. 10), and unpub- 
lished data by Weiland and Taylor were used in the present investigation. The average 
friction coefficient is calculated from the equation 


f_ 

2 


g^bav A Pfr 

4— G 2 

D 


(Bl) 


The friction pressure drop is obtained by measuring the total pressure drop and sub- 
tracting the momentum drop, 


Ap fr = Ap - Ap 


mom 


(B2) 


where 


Ap 


mom 


GRr2 _ h 

Sc\P2 Ply 


(B3) 


The average bulk density in the friction coefficient is evaluated at the average static pres- 
sure and the average static bulk temperature : 


o lfl +P 2 

tev ' 4 * 


(B4) 


Equation (B4) does not give the average density, which is 


p -L P JL + P J 

av Ah V 


(B5) 
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Equation (B4) is usually used to calculate the density that is used in equation (Bl). 
Equation (B2) can be written as 


A P = A Pmom + A Pfr < B6 > 

from which it can be seen that the accuracy of Ap fr is not only dependent on the accuracy 
of measuring Ap but on the magnitude of A P mom compared with Ap fp . Table II pre- 
sents the various investigators and the maximum ratio of momentum to friction pressure 

drop 4p mom /4 "’lr' 

The large ratios of AP mom to Ap fr (in ref. 10) resulted in an inaccurate Ap fr and 
correspondingly inaccurate friction coefficients. In some cases, the AP mom appeared 
larger than the total pressure drop across the test section, which resulted in negative 
friction coefficients. Because of these inaccuracies, the friction coefficients presented in 
reference 10 were not included in this investigation. 


Local Friction Coefficients 


The local friction coefficients presented by Perkins and Worsoe -Schmidt (ref. 4), 
Davenport (ref. 6), McEligot (ref. 7), and Magee (ref. 8) were calculated by the equation 



2g c 


(B7) 


where 


Tw = .5& + P vi 

4 dx 2g c 


(B8) 


All local friction coefficients seemed to have an acceptable degree of accuracy. 
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TABLE I. - TEST CONDITIONS FOR VARIOUS SOURCES OF DATA 


Investigator 

Ref- 

er- 

ence 

, Fluid 

Ratio of surface 
temperature 
to bulk 
temperature, 

V T b 

Friction 

coefficient 

type 

Length to 
diameter 
ratio, 

L/D or x/D 

Bulk 

temperature, T fo 

, Surface 

! temperature, T g 

Modified surface 
Reynolds number, 

Re s 

Tube 

orientation 

Tube inside 
diameter, D 


°R 

°K 

°R 

°K 

in. 


Taylor 

2 

Hg and He 

1.0 to 3.7 

Average 


829 to 


1559 to 

866 to 

1.9xl0 3 to 10.9X10 3 

Vertical 









1643 

913 

4749 

2638 





Taylor 

3 

H 2 and He 

1.0 to 4.1 

Average 



278 to 




Vertical 

.115 

.292 







1243 

691 


2555 





Perkins and 

4 

N 2 

1.4 to 7.4 

Local 

16,53,73, 


117 to 

371 to 



Vertical 

.124 

.315 

Worsoe -Schmidt 





93,113 




1127 





Davenport 

6 

N 2 and He 

1.1 to 2.2 

Local 

25,44,68 

721 to 




. 17X10 3 to . 94X10 3 

Vertical 

.125 

.318 







2226 

1236 


1366 





McEligot 

7 

He and Air 

1.4 to 2.5 

Local 

16,18, 27, 

541 to 


844 to 468 to 

5. lxl 0 3 to 40.9X10 3 

Vertical 

.123 

.313 






34,46,63 


568 

1705 

948 



.250 

.635 

Magee 

8 

He and Air 

1.4 to 2.5 

Local 

16,18, 27, 

541 to 


844 to; 468 to 

5. lxlO 3 to 40.9X10 3 

Vertical 

.123 

.313 






34,46,63 

1022 

568 


948 




.635 

Wolf 

9 

Air C0 2 

.35 to 2.7 

Average 

21,40,60 

545 to 

303 to 

544 to ! 302 to 


Horizontal 

.384 

.976 







1941 

1079 

1782 




.569 

1.445 

Weiland 

(a) 

Air 

1.0 to 1.5 

Average 


592 to 

329 to 

711 to 

395 to 


Horizontal 

.184 

.468 







1324 

736 

1896 

1048 





Taylor 

(a) 

He 

1.0 to 3.4 

Average 

80 


581 to 

3112 to 



Vertical 

.115 

.292 







1589 

883 








a Unpublished data, Lewis Research Center. 





TABLE n. - MAXIMUM RATIO OF MOMENTUM TO FRICTION 
PRESSURE DROP FOR EACH INVESTIGATOR 


Investigator iRefer- 


Fluid Maximum ratio 
I | of momentum to 

! friction pressure 
1 drop, 

Ap mom^Pfr 


Taylor 


Heating He 


McCarthy and Wolf 


Wei land 


Taylor 


Cooling Air 
Cooling C0 2 
Heating He 


Not given 
Not given 
.1 
.09 


^Unpublished data, Lewis Research Center. 
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